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We report a new possibility of using the McT2 (Constransverse mass) variable for mass measure- 
ment in single step decay chains involving missing particles with moderate transverse momentum. 
We show that its experimental feasibility is enhanced compared to the corresponding MT2-kink 
method. We apply this method to reconstruct a pair of chargino decay chains. 
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Many theories of new physics beyond the Standard 
Model (SM) are expected to provide a rich invisible en- 
ergy signal from their Lightest New Particles (LNP), 
which are stable Dark Matter candidates, missing in the 
detector of future hadron colliders. In this situation only 
the so-called 'M^ 2 -kink' method can provide the informa- 
tion necessary to determine the masses of both mother 
particle and missing LNP simultaneously for the events 
with a pair of single step decay chains (TJ [2] ■ 

The origin of the Mr2-kink is the variety of 'Extreme 
Kinematic Configurations (EKCs)' in the events which 
can contribute to the maximum of the Mt2 distribution. 
In general, for a given trial LNP mass, x, the different 
EKCs will take on different values, so in different regions 
of x, M™2*(x) will follow different functional forms. At 
the true value, x = m x, all of the EKC should provide 
the same Mt2 maximum value as the true mother parti- 
cle mass, my, by the definition of Mti 0- Consequently, 
M™ 2 ax (x) shows a slope discontinuity at the kink point, 
M- ax ( X = m x ) =m Y . 

Let us consider the system of a pair of single step decay 
chains at the LHC: p + p -> S T + Y 1 /Y 2 {-+ aX 1 //3X 2 ), 
where Y\ t2 are the two mother particles with identical 
masses, each decaying to visible a(/3) and missing LNP 
(X^)). St denotes the other remnants (and the sum of 
their transverse momenta) not from Yi decays, i.e. from 
the Initial State Radiation (ISR) or decays before Y{, so 
that they provide the total transverse momentum of the 
Y1+2 system of —St- In this event topology, there are two 
physical degrees of freedom that can generate the variety 
of EKC developing the kink. One is the invariant mass of 
visible particles m a ^ in iV(> 3)-body decays, generating 
a so-called 'Mass Kink' [1], and the other one is non-zero 
transverse momentum (St) leading to a 'Boosted Kink 
(BK)' 2 of the Y1+2 system of our interest. When a(/3) 
consists of a single visible particle and has a fixed invari- 
ant mass, then only the BK provides enough constraints 
for simultaneous mass measurement. 

However, the BK is not easy to identify [2] . It requires 
large St to have a clear kink structure. With the high- 
est practical values of St /my ~ 0(1 — 10) at the LHC, 
the kink structure is not clear enough to be reliably mea- 
sured. The profile of the St distribution is sharply de- 
creasing for large values at the LHC, and the endpoint 



structure of the Mti (x) distribution is smeared for large 
St with small statistics. This can introduce significant 
systematic errors in the endpoint extractions. 

In this paper we present the endpoint behavior of the 
Mct2(x) distribution [4] for single step decay chains with 
non-zero St- The Mq^(x) distributions for different St 
carry independent information about my and mx, and 
it can be utilized for simultaneous mass measurement of 
mother and missing particles. Furthermore, as the de- 
pendence of McT2 on St is much more significant com- 
pared with the Mt2, it becomes good collider basis in 
resolving the boundary of the distorted phase space by 
various transverse boost effects. One more interesting 
property is that the endpoint structure is statistically 
cleaner with much less pollution from unwanted width 
effect in addition to the amplification and sharpening of 
the structure as noticed in [3]. In next sections, we de- 
scribe the merit of using M™|^ over M^g* for systems in 
which reasonable St is expected from heavy colored NP 
particle decays. We demonstrate that it can be used to 
probe the boosted kinetic boundary, and determine the 
masses of the particles involved in the chargino decays 
into sleptons, in particular X\ —> ^(^ Xi + v ) + 

The Mct2 (Constransverse mass) variable [4] for the 
Y\ +2 system is defined by replacing M T with M CT [5] in 
the definition of Mt2 [3] as follows: 



M C T2(X) = min 

klT + k2T — p T 

r(l) 



max {M^,M^} 



L CT> 1V± CT 

M ( c T (xY = x" + 2(|a T | ei +a T -fciT) 
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where Mq^,(x) is defined by replacing (ay, ex, kyr) with 
(hi G2i fer)- fci,2T and ei t 2 are the transverse (Tr) mo- 
menta and Tr energy, respectively, of missing X\^ 2 , with 
total missing Tr momentum, Tf, T = — («t + Pt) — St- 
o-t(Pt) represents the Tr momenta of visible particles 
from 11(12) decays, — St is the total Tr momentum of 
the Yi+2 system, and x denotes the trial test mass of Xi. 
Visible particle masses are assumed to be zero, however 
the definition of Mqt2 can be easily generalized for mas- 
sive visible particle events. In 0], it was found that the 
number of events near can be amplified with a judi- 

cious choice of trial mass, x, and Mct2 was employed to 
measure some physical constraints involving squark and 
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gluino masses with fewer systematic errors in pinpoint- 
ing the endpoint than with the Mti analysis. One thing 
worthy to note here is that interestingly, the endpoint 
of the Mti (x) [H 12] or Mqt2 (x) 01 f° r anv X values, is 
invariant under the unknown \fs{> 2my) values of the 
system with two mother particles, and this property is 
also satisfied by the M^^ix) f° r general boosted decay 
system as in following paragraphs. 

When St = 0, the Mct2{x) nas a one-to-one corre- 
spondence with Mti (x) j an d its maximum value is given 
as follows H] : 



M^f 



(X) =X 2 + 2(ao£f 



(2) 



where ao 



and E° x 



Vx 2 



ai. However, 



if St ^ 0, simple one-to-one correspondence no longer 
exists, except for the events of some EKCs. Certainly, 
depending on %, Mct2{x) an d Mt2(x) can be one-to- 
one, however, it's not valid in general x values. In result, 
the maximum profile of the Mct2(x) distribution shows 
a 2nd order Boosted Kink, ( 2 BK), structure which is dif- 
ferent from the BK of Mt2 ■ As investigated in [TJ [2] , the 
M™2 X (X] $t) value of the Y1+2 system is the same as the 
M^ 0:c (x,<5t/2) of single F* decay system, and the EKC 
for the M™2 X can be characterized by a pair of identical 
EKCs corresponding to the M™ ax of single Y^ decays al- 
though they have to be combined to produce the event 
with general y/s > 2my . Then, in the language of single 
Yi decay, the BK of M^ x {x) is provided by two EKCs 
of the single Y t decay events, characterized as follows : 



for x < m x b) 



7T for x > m x 



where (\> max is the azimuthal angle between visible Tr mo- 
mentum and St in the rest frame of the Yi with vanishing 
Ar)(= r) vis -r) inv ). Similarly, the EKCs for Mg^f (x) with 
2 BK can be also characterized by single Yi decays, but 
the azimuthal angle dependence is different : 



1) COS( 

2) cosi 



\S T \E Y (m Y -m 2 x ) 

-1 for x > X*, 



+ 



\s T \ 2 



for x < X* 



where |<5 T | = \8 T \/2 and E Y = \jm\ + |<5 T | 2 . Thus, 
McT2.(x) is described by two curves in the two regions 
of x divided at x*> where x* is given by 



xl = \~6 T \{2a-\S T \), 



(3) 



with a = |a |(|r + 1^)- Then M CT2 for thc ful1 X 
range is given as follows : 



M^f 



for x < X* 



(4) 



2 X 2 + \St\ 2 

= x 2 + M\St\ 2a^ X 2 + (|*r| - «) 2 

forx>X*. (5) 

When x* is real, the two maximal curves of Eq. Q and 
([5]) come into contact at x* with the same inclination, 
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FIG. 1: a) M^rUx) for (m Y ,m x ) = (130, 100) GeV and 
St = 247.9 GeV with 8 T = 20 GeV (x* = 21.8 GeV) b) 
M C T 2 (x) distribution for 8 T = 250 GeV (No 2 BK). 

but the curvature of M^^ix) nas a discontinuity at x = 
X* which we denote 2 BK. It's interesting that when a 
2 BK exists, there is a boost-trapped distribution for x < 
X* in that its boundary is independent of any physical 
masses in the decay system, as shown in Eq. Q. If x* 
is imaginary, Mgjf(x > 0) just follows the line given by 
Eq. ([5]). For given parameters, (mY,mx, \8t\), X* is rea l 
and observable if 

1. < ^, then x* G R. 

2. < ^ < 1, then X^GR if \S T \<\8 T \, 



V z "by 

where \8 T \ = 4a / ' ^Jl~ 



4a 
my 



The reality condition of x* for the existence of a 2 BK 
is not always met, however, it certainly appears for the 
events with small values of \St\- For example, Fig. ([I]) 
shows contour histogramming on the Mct2{x) ~ X vs X 
plane. We used the 10 5 Monte-Carlo events which simu- 
late the phase space of the pair of boosted decay system, 
Y 1+2 + 8 T where (m Y ,m x ,8 T ) = (130,100,20) GeV with 
massless visible particle. In Fig. Qa)) the M£?|^(x) 
is well described by the two curves in the two regions 
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of x divided by 2 BK, x* = 21.8 GeV. For X < X*, 
Mcriix) follows the blue curve; it switches to the red 
one for x > x* with a continuous slope at x — X* ■ Since 
the spectrum for Fig. ([I]) belongs to the second cate- 
gory of reality condition with S T — 247.9 GeV, no 2 BK 
arises for St — 250 GeV > S T as shown in Fig. (Urn)). 
In any case one can try to fit the endpoints with the 
segmented maximal endpoint function, Eq. Q and ([5]), 
for given maximal value of \St\, \S T lax \ in an event sam- 
ple, since the M™|^(x) always increases with \St\- The 
segmentation point, x* °f the piecewise fitting function, 
Eq. Q and ^ is also a function of a, and overall fit- 
ting procedure systematically fits on a without a prior 
knowledge about x*- Then for several |#p ax | values and 
given x values, one can extract independent mass con- 
straints, a(my, mx, |#p ax |) by which the my and mx 
can be simultaneously measured. 

The Mct2 is found to be so sensitive under the Tr. 
boost of the system that even with moderate values of 
non-zero <5t, it can provide good resolving power in dis- 
crimination of various transverse boost effects which are 
usually buried in the large systematic errors of shifted 
endpoint extraction. The rates of M™£y T2 (x) shift with 
respect to a change in St are given as follows : 



F) M max 

U1V1 CT2R/T2R 



d\s 2 



(6) 



aE^ 



f max 

1 CT2R/T2R* 



1 ± 



1*1 



1± 



for M CT 2/T2i respectively, with E x = \/x 2 + (|*r| ~ a ) 2 
and Mt2R = X 2 - 2a(|*r| - a) + 2aE x . Here the 
subscript letter- R denotes the values of M CT2 /T2(x) f° r 
the x values larger than corresponding kink positions, 
X* / mx ■ The shift of the is magnified in compared 

with M™ as indicated from the flipped signs between 
the rates of the endpoint shifts in Eq. In particular, 
the shift can be large by the factor of O(IO) 1 , when the 
mass difference between My and Mx is small enough so 
that \&r\ > ol for a moderate value of \S~t\- in Fig. (|2la)) 
we show the shift of the M™|^ compared with that of 
M™ 2 ax under the effect of non-zero St- We define R as 
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FIG. 2: a) Ratios of Sr-shift between Mg^l and M£f x with 
respect to St /My for various mass spectrums, Mx/My- b) 
<5 T -shift of Mg 1 !! - M^l{m x = 0) for m x = (10 - 190) GeV 
with the do = 62.4 GeV and the trial LNP mass, x = X* + 
50 GeV. 
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We take A = my/2. Also, if mx/my = 0.1 
Xi = 2x*, where x* is calculated with \St\ = |*r 



0.7, 

-A, 

while xi — m x + my/2 if mx/my — 0.8 — 0.9. Un- 
der the transverse boost effect, the magnified shift of 
^ctHx^t) can be clearly seen depending on the var- 
ious mass ratios of mx/my in Fig. (121 a)). If we com- 



1 Actually, the sensitivity is not only for the EKC events. The 
relation |6j| holds for all events with general visible transverse 
momentum smaller than a. 



pare the shifts of the ^ct 2 /T2 using specific mass spec- 
trums with transverse boost effect, (my,mx, |*r|, A) = 
(150,100,100,100)|(500,100,100,100), then the shifts 
of Mg iax /T2 (x) are calculated to be AM™ ax (x = 
X*) = 73.5|38.4GeV while AM™ ax (x = 200 GeV) = 
12.7|11.6GeV. 

In calculating the R, the optimal choices of xi,2 for 
which each of the shifts is maximized and well-measured, 
are ambiguous. We choose xi ~ 0(x*) because it al- 
ways provides sizable shift with clean and sharp end- 
point structures of the Mqt2 distributions [8]. A prac- 
tical choice of X2 is also unclear. Indeed, as we take 
larger value of X2j the denominator of R increases, re- 
ducing R down to 0(1) in some cases. However, as we 
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will see in next section, there appear faint long tails near 
the expected endpoints of the Mti distributions for large 
St, making it difficult to measure the endpoints and the 
shifts. Also, for a given value of St, such an endpoint 
structure becomes worse as the x increases over the true 
value of x = m x , while the denominator of R gets rapidly 
saturated just with x ~ 0(mx) > mx- Thus for the 
sake of both of the endpoint sharpness and the sizable 
endpoint shift, we choose X2 ~ m x ~ my/2 so as to fix 
the separation of x from the true mass, mx to be the 
order of mx, at least to be m x /2 in Fig. (j2ja)). 

Fig. (||») displays the <5 T -shift of Mgffi - 
M™^2( TO x = 0) f° r various missing LNP masses, mx = 
(10 - 190) GeV while the a is fixed by 62.4 GeV and 
the trial LNP mass, x = X* + 50 GeV. For St — 0, no 
mass resolving power is obtained by the endpoint mea- 
surement as expected, however, one can see that the St~ 
shift of M™|^ can be large enough to measure the new 
particle masses with the resolution of 0(1 — 10) GeV for 
St ~ O(10 — 100) GeV. It is also enhanced for the large 
mx I my case. We fixed the ao because it is the basic 
momentum scale one can observe in detector regardless 
of the mass spectra. 

By projecting events in the Mct%{x ~ 0(x*)i $t ^ 0) 
basis, one can get a more sharper edged and (J^-sensitive 
event distribution. The flipped sign in the definition of 
the Mcti variable makes the distribution compact with 
respect to the internal momentum scale of the system, 
while being much more sensitive for the external boost 
momentum St, like a rubber ball. This means that one 
can have a better chance to measure both of the masses 
with enhanced resolution in probing the boosted kine- 
matic boundaries. 

In the following example, we employ the Mcti (x) vari- 
able for resolving superparticle masses, especially in the 
events with a pair of Xi decays. It is well-known that 
one can efficiently select such NP events by imposing the 
same-sign (SS) dileptonic condition 9J. This signature 
provides an important discovery channel in establishing 
SUSY, because the background is small, however, mass 
measurement in this channel is difficult because of the 
existence of extra missing vs and combinatorial SUSY 
backgrounds of Xi (—> + + + decays. There 
exist a few studies [5] on mass measurement in this event 
topology. We demonstrate that the mass resolving power 
is significantly enhanced by using the M™ T ^(x) under re- 
alistic \5t\ values of Xi pahs arising from squark/gluino 
decays. 

SUSY example We demonstrate the mass measure- 
ment of X\ and v for a SUSY benchmark point de- 
scribed in Ref. [7J. In this scenario, 3rd generation 
sfermions are very heavy and X\ has no branching ra- 
tio (BR) to t with BR(xi -> I + v) = 0.63 fol- 
lowed by v — > Xi + v ■ The relevant mass spectrum 
of the point is (m^m^m-i^^m^) = (~ 720,^ 

635, 231.5, 157.2, 123.3) GeV with a^ SY ~ 17 pb. We 
calculate Mcti (x) f° r a system of a pair of SS chargino 




FIG. 3: Contour plots of the SS SUSY dileptonic events sur- 
viving the event selection cuts in Mqt2{x) ~ X vs X plane for 
|5£ lax | = (50 T op-Lcft, 160 T -Ri g ht, 300 B otto m -L, 400 B - R )GeV. 

decays, where ctT and Pt are the transverse momenta of 
leptons from the chargino decays, and the two missing 
/>s' momenta provide If, T [10] . All of the other jets in the 
event are taken as the St of the chargino decay system. 
SUSY events of 70/6" 1 at ^fs = 14 TeV are generated 
by PYTHIA [11] for this point. The detector simulation is 
implemented with PGS4 [TJ] . 

We select SS dileptonic events with N(> 2)-jets + 
Jp T . The event selection cuts are Pt{£i,2) > 20 GeV, 
Pr(jet h2 ) > (100, 80) GeV, and # T > 100 GeV. Af- 
ter the cuts, the signal to background ratio becomes 
high and the kinematic endpoints are expected to ap- 
pear clearly without background from the SM events 
and neutralino decays [5J. Additionally, we impose a 
cut on the \S T \ value, |&r|< |#? ax |- For S iven l#r ax | 
and x values, M™|^(x) ( 4j5 | of the signal events is 
completely described by two known (\St\ — \S- 
and two unknown (m~±,mv) parameters. Fig 
shows the contour histogram of the events which sur- 
vived the cuts in the Mct2{x) ~ X vs X plane for 
|#p ax | = (50 (a) , 180 (b) , 300 (c) , 400 (d) )GeV. The events 
are p opulated below the expected endpoints of Mct2 (x) 
( |4j5] ). The solid lines show M£?||(x) @ for all of 
the |(5™ ax | values. The boundaries of the Mct2 dis- 
tributions are well described by the solid lines in the 
range of X < X*{= 53.4 (a) , 91.3 (b) , 109.3 (c) , 120.7 (d) ) GeV 
in Fig. Q. In each plot, we also show three func- 
tional lines of Mg™( X > X*,\5t^\) © for three dif- 
ferent mass spectrums, TO y(DashodiDashod _ Dottcd!Dottcd) = 
(m^± , m - ± — 30, m~± — 60) GeV, with the same mx = m„ 

X± Ai Xi 

values. The boundary of the Mqt2(x, I^T ax l) distribu- 
tions also appear to be well-matched with the expected 
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FIG. 4: a) T " L |<5 T | distribution of the SS dilepton + $ T 
system surviving the cuts. b) T ~ R Mct2 and c) B " L Mt2 distri- 
butions. d) B " R Consistent mass ranges in My — Mx vs Mx 
plane. 

Dashed-lines for x > X*- 

Using the reconstructed Mq^(x) Ior different |£™ ax | 
values, we finally attempt to get the consistent range of 
chargino and sneutrino masses. Fig. (Ha)) shows the 
\St\ distribution of the SS dilepton + Ipj, system after 
the cuts. The event statistics near \6t\ ~ O(100) GeV 
are very rich due to the large boost effect of hard jets 
from initially produced colored superparticles, and it can 
be a good testbed for the mass measurement utilizing 
various transverse boost effects. For different values of 
\6™*\ = (400 S oiid,300Da S hcd,200Dottcd,50 D -D)GeV cuts, 
corresponding Mcriix = X* + ~ 20 GeV) distributions 
are shown as in Fig. Qb)). Fig. Qc)) displays 
M T2 (x = 400 GeV) distributions with the |<5f? ax | cuts. 
In both of the plots, red arrows denote the expected end- 
points of the M C T2/T2{x) f° r the given |<5™ ax | cuts. Then 
one can clearly see that 

1. The endpoint structure of the Mct2{x) distribu- 
tion is statistically cleaner with less pollution [8] 
and sharper [3] in wide range of |<5™ ax |, compared 
with Mt%{x)- 

2. The Mqt2 is so sensitive under the transverse 
boost effect, the shift of M™|^ with a change of 



|<5™ ax | can be sizable well beyond the uncertainties 
in endpoint extraction, as well as the number of 
events contributing to M CT2 (x,\^t\ > \S'^ ax \) > 
^CT2(X: I^T ax l)' i s mucn larger than the num- 
ber of the events with M T2 (x,\ s t\ > \5!p ax \) > 

M?p(x,\sr x \) 

Putting all these facts together, we now have a collider 
variable with which the measurement of Tr. boost-shifted 
endpoints are experimentally reliable for mass determi- 
nation under various systematic uncertainties of endpoint 
extraction. 

An error profile function of x, — a \ £>(x — X*) + c 
is assumed in the two regions of (x < X* \ X > X*)j 
which are likely to be obtained by simple endpoint fit- 
ting of Mct2{Xi I^T ax l) distributions. The coefficients 
are conservatively set to (a, b, c) — (2|3, 0.03, 2|3) GeV for 
|#p ax | = (50|400) GeV. Then, for the two |#p ax | values 
we fit two a values in Eq. ([3]) and ([5]) , using the parame- 
terized endpoint measurements, M™J^(x, |<5™ ax |)=l=e(x) i n 
the range of < x < 300 GeV, including \* which is also 
a function of a. This is a one-parameter fitting of a for 
a given |<5™ ax | value, and the precision of fit result is not 
much affected by a change of the fit region of x, as long as 
it is taken broadly (0,x* +O(50— 100) GeV) enough to in- 
clude the x* region. For the parameterized errors, a±5a 
values are fitted to be (69.5| 136.5) ± (3.38|5.14) GeV for 
\S™ ax \ = (50|400) GeV, with x 2 /ndf < 1. Fig. gd)) 
shows two bands(Greon— Black) in My — Mx vs Mx plane, 
consistent with the fitted a values. The black dot indi- 
cates the true mass point and the yellow intersection re- 
gion is the lcr range of the mass measurement using the 
M CT2(x) data for two different \S^ ax \ cuts. 

Discussion In this paper we study the Mct2(x) end- 
point of a system containing a pair of single step de- 
cays, which is expected to be the simplest non-trivial 
new physics event topology. The Mqt2 magnifies the 
transverse boost effect on the system, while keeping 
its endpoint structure clean and sharp. Thus, measur- 
ing Mqt2 f° r different \St\ values can provide us with 
reliable resolving power to determine both masses of 
mother and missing daughter particles only using mod- 
erate transverse boost momentum of the system of our 
interest. With realistic St profiles of mSUGRA like events, 
we demonstrate simultaneous mass measurement of Xi 
and v in SS dileptonic decay channel, where xt v + 1. 
The technique also works for a single mother particle de- 
cay system with non-zero transverse momentum. 

This work was supported by World Premier Interna- 
tional Center Initiative (WPI Program), MEXT, Japan. 
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